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"ob"; _p 

. We investigate the — D mixing through the doubly Cabibbo suppressed (DCS) 

^ ; channel D° /o(980)i^*° and its charge conjugate channel, in which the K 

^Cp ' meson is reconstructed in both K^tt and Ksir^ final state. Although the decay 

QQ ' /o(980)i^* has a small branching ratio, the final state mesons are relatively 

easy to identify. The /o(980) meson can be replaced by tt'^it~ in which vr+vr^ form 

' an S'-wave state, or a longitudinally polarized vector meson p^,uj. All mixing pa- 
I ' 

pi i' rameters, including the mass difference and decay width difference, can be extracted 

. by studying the time-dependent decay width of these channels. We show that the 

I method is valid in all regions for mixing parameters and it does not depend on the 

> ■ 

i/^ ' strong phase difference. 

^ : 

in ■ 

^ ' The neutral meson mixes with its CP-conjugate D° through box diagrams in the 
CN ; 

O I standard model (SM). In the box diagrams, the b quark coupling is small and the s and d 

^ ' quark have small masses. Thus the mixing oiD^-D in the SM is very small, measurements 



X 
J3 



on the mixing parameters serve as a probe to detect new physics scenarios. In the presence of 
mixtures, the interference in transition amplitudes of D and _D° is possible. Experimentahsts 
can study the time-dependent decay widths to extract the mixing parameters. Recently, 
e and CDF collaborations have reported their measurements on — 



the BaBar, Bel 
mixing 



QlflflB 



5| . Meanwhile, a number o: 



Id, 



11 



theoretical methods are proposed to measure 



13|. 



the mixing parameters 

In this note, we propose a new method to extract the mixing parameters in — 
mixing: using the doubly Cabibbo suppressed (DCS) channel —>■ MK*^ (M denotes a 
/o(980) meson, a non- resonant S'-wave vr^vr" state, or a longitudinally polarized p° or u 
meson) and its CP-conjugate channel. The decay MK*^ can proceed in term of 

the DCS amplitude or the amplitude from the mixing followed by a Cabibbo allowed decay 
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MK*^. The two amplitudes interfere with each other. If the K*'^{K*^) meson is 
reconstructed in both K^tt^ (K^tt^) and Ksir^ final state, the mixing parameters can be 
determined through the time- dependent studies on the decay widths. 

At present, there is no direct experimental measurement on —>■ fo{980)K*^. The 



branching ratio for —>■ fQ{980)Ks is given by 14l |: 



BTZiD^ ^ fo{980)Ks) = {IMt'oH) x 10-^ (1) 

where /o(980) is identified in the vr+vr" final state. This value can be used to estimate the 
decay rate of — > /o(980)/i*°, together with the following results |l4| : 

Bn{D^ ^ Tf^Ks) = (1.14 ± 0.12)%, 
S7^(Z}° ^ 7r°ir^°_^+) = (1.91 ± 0.24)%, 



where the subscript denotes the daughter mesons to reconstruct the K*'^ meson. Under the 
factorization assumption, the decay amplitudes are proportional to ^°(/o) form factor 

times the Kq{K*'^) decay constants. From Eqs. and ([2]), we expect a somewhat larger 
branching fraction for the D° fo{980)K*^ decay channel than Bn{D° fo{980)Ks). 
Since the sum of masses of /o(980) and K*^ is close to the mass of meson, this decay 
channel will be suppressed by the phase space. To obtain some information on the phase 



space suppression, we can compare the two decays D° — > (j)Ks and D° (pK 



i37^(D° (pK+K-Ks) = (2.10 ± 0.16) X 10-^ 
Bn{D° (f)K+K-K*j^~^+) = (1-01 ± 0.20) X 10"^ (3) 

The branching ratio of —>■ (j)K*^ is about 5% of that of —>■ (f)Ks- Assume the same 
suppression for fo{980)K*^, we expect this decay has a branching ratio of O{10~^) 

which can be studied by the future experiments. In our method, one can also replace the 
/o(980) by the non- resonant tt+tt^ final state where vr+vr" is restricted to be an S-wave state. 
The three-body (non-resonant) decay —>■ 7r+7r^/i*° also possesses a sizable branching 
fraction [14]: 

Bn{D° 7r+7r-i?^°_^+) = (9.7 ± 2.1) x 10^1 (4) 
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In the following, we will take D — > /o(980)ir* as an example and use /o to denote /o(980) 
for convenience. 

For foK*'^{K*'^), the transition in quark level is either c sdu or c ^ dsu. The 

former transition is proportional to VcsV*^ ~ 1 and the latter is suppressed by the CKM 
matrix elements: VcdV*^ ~ 0.04. Since there is only one amplitude for each decay in the SM, 
the direct CP asymmetries are and the amplitudes satisfy the relation: A{D^ /qK*^) — 
A(D° foK*°) and A{D^ foK*°) = A(D° foK*°) = Af^K*- Neglecting direct CP 
asymmetries, one often defines the two parameters rf^K* and Sf^K* by: 

Under the assumption of CPT invariance, mass eigenstates of neutral D meson system 
are given by: 

\Di) = pm + ql^), \D,) = p|L>0) -q\^), (6) 

where p and q satisfy the normalization condition + = 0is defined as the phase 
of q/p: (f) — arg(g/p). In the SM, the phase is very small and thus the large value for 
this parameter will definitely imply the presence of new physics. Di and D2 have different 

masses and different decay widths. The differences are defined as: AM = Mi — M2 and 



AM 

r ' 



AP = Pi — P2. For convenience, experimentalists also use the two parameters x 
y = where P is the averaged decay width of the L>- mesons. 

In the limit of x <^ 1, y <S 1 and Ft <^ 1, the time-dependent decay amplitudes squared 
of D° — > / in which there is a purely D° state at i = is given by: 

\A{D\t) - /)r = e-r* [Xf + YfTt + Zf{Ttf + ...] , (7) 

where the expansion is only at 0{TtY accuracy. It is similar for \A{^{t) /)|^: 

\A{-^{t) = e-r* \Xf + YfVt + Zf{Ttf + ...] . (8) 

In the following, we will take / = f^K*^ and / = f^K*^. If the K*^{K*^) meson is recon- 
structed in K^n^ mode, the coefficients in the time-dependent decay amplitudes squared 
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are given by: 



^foK* = \-\\AfoK*\'^rf^K*{y'f,K* COS0 - x'j^K. sin0), (10) 
YfoK* = \-\\Af^K*\'^r f^K* (y'foK* COS (p + x'f^K* sin (p), (11) 

ZfoK* = |-| |^/oK*| ^ , (12) 

ZfoK' = \-\ \AfoK'\ — ^ — • (13) 
x' and are linear combinations of x and y: 

x'foK' = i.x cos Sf^K* + 2/ sin 5foK* ) , 

?//oi^* = cos - a; sin 5f^K* )■ (14) 

The amplitude squared |y4/QX*P can be easily measured using the time-integrated rate for 
Cabibbo favored mode foK . In foK , the Y and Z terms come from the 

mixing followed by the DCS decay D f^K which are at least suppressed by Vf^^K* ~ 
0.04. Neglecting the power suppressed terms, only the first term X = \AfQK*\'^ survives and 
the time-integrated rate is given by: 

\A{D'{t) - foK*')\'dt ^ (15) 

With P determined, the ratio rj^K* is easily determined by combing the measurements 

on XfgK* and |y4jgX*P- Furthermore, \q/p\ and x^ + y"^ are determined using Zj^k* and Z j^k* 
in Eqs. f fT2|) and f|T3|) where terms proportional to rj^x* neglected. 

Using the decay widths given in Eq. ([7]), Eq. ([8]) and Eq. (fT5l) . the four parameters are 
determined while the other ones (0, Sf^^x*, x/y) still remain unknown. The decay — »• 
foK*^ where the K*^ is identified by the Kstc^ final state could provide more measurements 
which are helpful to extract the other parameters. The amplitudes receive two parts of 
contributions: 

A{D' ^ /oi^x,.o) = A{D' ^ UK*x,^o) = Aj,kA^ - rf,K*e-''w^*] ^ Af.Ks.. (16) 

where the appropriate reconstruction factors with a K*^ identified by the Kst^^ are assumed 
in the experimental measurements. The coefficients in time-dependent decay widths are 
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given by: 

XfoKsTT = XfoKs-K = \AfoKsn\^, (17) 

YfoKsTT = -\-\\Af^,Ksnfi-xsm(j) + ycos(j)), (18) 
^foKsn = -\-\\Af^^Ksn\^{xsm(f) + ycos(f)), (19) 

where we have added the subscript KsT!' for the coefficients. With these measurements, the 
mixing parameters, tan^ and (x/y)^, are solved by: 



tan = = — , (20) 



o T. 

X 



(21) 



where cos2(/) is chosen positive (very small (j) in SM). For convenience, the two coefficients 
have been reexpressed by: 

y^L = ' " :i^!r""' = (22) 
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2X,,,,.|g|b| - ^'^^ 

and J^foKsn = P + Yft^Ks-K ^ ^/oi^sTT- the limit of ^ 0, the measured value of Y^^j^^^ is 
very small and thus we can expand the above solutions in Eq. (l2Up and Eq. (12 ip into: 

tanV= -^^ + ..., (24) 

2 /•2_y(+)2 y(-)2 .2 

y2 ~ ^(+)2 ^(+)2 /-■2_y^(+)2 ^^^^ 



where higher powers of Y^^j^^^ were neglected. As we can see from Eq. (!25|) . there are two 
advantages in our method: the ratio of and y"^ is finite even for tiny 0; our method does 
not need the strong phase difference 6 which is unknown at all. 

The authors in Ref. [3] propose to extract the mixing parameters using the time- 
dependent study of K*'^TT^. The four parameters, [^^^(^o)^^. p, \Af^(^T,o)K*\'^, \p/q\ 
and p = + y'^, are determined in the similar way but the other parameters are extracted 
in different ways. In —>■ Kstt^tt^ decays, the final state contains two neutral pion. It 
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requires four photons to reconstruct them which is a rather difficult job. Thus they suggest 
to use the normahzation of Kst^^t^^ to give a constraint on the strong phase difference. 

In the present method, the normahzation constraint is not used and we utihze measurements 
on the time-dependent decay width of foK^^^o instead. Ahhough the foK* 

has a smaU branching ratio, the advantage is that the /o meson is easier to re-construct than 
7r° on the experimental side. 

In summary, we have studied the D'^ — mixing in D ^ foK* decay, where the K*^ 
meson is reconstructed in both K^7t~ and Kstt^ final state. The method can be directly 
generalized to the three-body decay D tt^tt^K* where the vt^tt^ is an S-wave state. The 
/o meson can also be replaced by a longitudinally polarized p° or u meson. All parame- 
ters in D'^ — mixing can be extracted by studying the time- dependent decay widths of 
these channel. The extraction is valid in all regions for mixing parameters and it does not 
depend on the strong phase difference defined in Eq. ([5]). The future experimental stud- 
ies on D ^ foK* decay, together with the D p^{uj)K* decays in which vector mesons 
are longitudinally polarized and the three-body decay D n^n^K*, can provide another 
alternative method to extract the mixing parameters. 
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